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We study scaling behavior in 2D, square lattice, S = 1/2 and S = 1 Heisenberg antiferromagnets using the data on 
full q-dependences of the equal time structure factor and the static susceptibility, calculated through high temperature 
series expansions. We also carry out comparisons with a model of two coupled S=:l/2 Heisenberg planes with the 
interlayer exchange coupling tuned to the T = critical point (two-plane model hereafter). For both S=l/2 and S=l 
models, we separately determine the spin-wave velocity c and mass ni = c/^, in addition to the correlation length, ^, 
and find that c is temperature dependent; only for temperatures below T ^ JS, where J is the exchange coupling, c 
approaches its known T = value, cq. This non- universal lattice effect is caused by the quantum nature of spin, and 
is therefore not captured by the quantum nonlinear a-model. Despite this temperature dependence of the spin-wave 
velocity, full q~ and oj-dependences of the dynamical susceptibility x{^.^ ^) agree with those of the universal scaling 
function, computable for the cr-model, for temperatures up to Tq ^ 0.6co/a, and their further analysis leads us to 
the inference that below Tq the S=l model is in the renormalized classical (RC) regime, the two-plane model is in 
the quantum critical (QC) regime in agreement with earlier work, and the S=l/2 model exhibits a RC-QC crossover, 
centered around T=0.55J. In particular, for the S=l/2 model above the RC-QC crossover and for the two-plane model 
at all temperatures where calculated, the obtained spin wave mass m = c/S, is in excellent agreement with the known 
universal QC prediction, m ~ 1.04 T. In contrast, for the S=l/2 model below the RC-QC crossover, and for the 
S=l model at all temperatures, the behavior agrees with the exact RC expression. For all three models, nonuniversal 
behavior occurs above T 0.6co/a. Our results strongly support the conjecture by Chubukov and Sachdev that the 
5 = 1/2 model is close enough to the T = critical point to exhibit QC behavior at intermediate temperatures. 
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I. INTRODUCTION 

The problem of antiferromagnetism in two dimen- 
sions has recently attracted much attention, in part 
due to its relevance to high temperature supercon- 
ductivity. Since the high temperature superconduct- 
ing materials are derivatives, obtained by doping, of 
parent insulating layered S=l/2 antiferromagnets, the- 
ory of the insulating phase is important for under- 
standing the unusual magnetic and transport proper- 
ties in the normal state of doped materials. Our re- 
sults are directly relevant to the following two materi- 
als: La2Cu04, which is the insulating parent compound 
of the La2-xSrxCu04 and La2-xBaxCu04 superconduc- 
tor families, and Sr2Cu02Cl2, which is unrelated to any 
superconducting family, but has well studied magnetic 
properties, similar to those of La2Cu04. Both materi- 
als experience antiferromagnetic transitions to the Neel 
state with nonzero staggered magnetization, La2Cu04 at 
Tn =i 300K and Sr2Cu02Cl2 at Tn ^ 256K. According 
to the Hohenberg-Mermin- Wagner theorem, in genuine 
two dimensional antiferromagnets the long range order at 
finite temperatures is prevented by thermal fluctuations. 
In real materials, small interlayer exchange interaction al- 
lows the transition to occur at a finite temperature Tjv. 
In what follows we discuss only the properties well above 
T/v, where the magnetic subsystem can be considered 
two-dimensional. 

It is well established [|] that the magnetic properties of 
La2Cu04 and Sr2Cu02Cl2 are described by the S* = 1/2 
Heisenberg model on the two-dimensional square lattice. 
The Heisenberg model is defined by the following Hamil- 
tonian: 



H 



(1) 



where S^ are local spin operators and the notation 
means that the sum is taken over all pairs of nearest 
neighbors. In what follows we set ks = h ~ a ^ 1, where 
a is the lattice constant, although in certain expressions 
we retain a for clarity. The dynamical spin susceptibility 
x(q, a;) = x' + ix", is defined such that the magnetiza- 
tion induced by a wavevector and frequency dependent 
magnetic field is: 



M(q,w) = x(q,'^)Hcxt(q,t^) 



(2) 



We calculate high temperature series expansions in pow- 
ers of /3 = J/T for two quantities: the static q-dependent 
susceptibility, defined as a response to a static, but gen- 
erally nonuniform, magnetic field. 



x(q) = x(q,^^ = 0), 



(3) 



and the Fourier transform of the equal time correlation 
function of two spin operators. 



5(q) -^(S,(0)S,(r))exp(iqr), 



(4) 



where angular brackets represent thermal averages. 
These quantities can be expressed in terms of x"{l,^), 
which plays the role of the spectral weight for spin fluc- 
tuations, according to: 



n— — OO ^ 

1 

x(q) = x(q, iujn = o) = - c 

J-oo 



th(w/2T)' 
X"(q,^) 



(5) 



where LUn = 2TmT is the imaginary bosonic Matsubara 
frequency. 

Total spin conservation in the model leads to the rela- 
tion 



Xo = ffVlx(q = 0) = g'^ilT-'S{q = 0), 



(6) 



where Xo is the bulk magnetic susceptibility, g the elec- 
tron gyromagnetic ratio, and the Bohr magneton. 

We study both the case of S" = 1/2, which corresponds 
to La2Cu04 and Sr2Cu02Cl2, and the case of = 1. 
We also study a model of two S=l/2 coupled Heisenberg 
planes at its critical point, J2/J1 — 2.51 ± 0.02 (2|. For 
both x(q) and S{q), we are able to generate series for 
an arbitrary wavevector q complete to order for the 

= 1/2 model, and to for the S" = 1 model. For 
the two-plane model we set J2/J1 = 2.5, and generate 
series for the structure factor and the susceptibility at the 
antiferromagnetic ordering vector to order Z?^". We use 
Fade approximation techniques for series analysis. We 
emphasize that the series converges extremely well in the 
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temperature region presented in the paper. In fact, with 
the available long series, most of the Pade approximants 
are within a per cent, or so, from each other, down to 
T « 0.4 J for 5' = 1/2 and T « 1.5 J for S = 1. This 
allows us to put quite small errorbars on our results for all 
but the lowest temperatures. At the lowest temperatures, 
the approximants begin to significantly differ from each 
other and hence we do not go to still lower temperatures. 



II. OVERVIEW OF RENORMALIZED 
CLASSICAL VERSUS QUANTUM CRITICAL 
SCALING BEHAVIOR 

Scaling behavior of the 2D Heisenberg models follows 
from their mapping |^ to the quantum nonlinear sigma 
model (cr- model hereafter). When the long range order 
is present at T = 0, fluctuations at low enough tempera- 
tures are purely classical |^ , in other words the character- 
istic frequency of spin fluctuations uj <^T . This happens 
because for T ^ ps {ps is the T=0 spin stiffness) , 



T^/^e^^{--2nps/T) 



(7) 



decreases faster than the temperature as the tempera- 
ture decreases. Accordingly, this regime is called the 
renormalized classical (RC) regime, where "renormal- 
ized" means that the T ~ parameters, such as the 
spin stiffness, ps, and spin wave velocity, c, are renormal- 
ized compared to their mean-field values as a result of 
quantum fluctuations at shorter wavelengths. 

Recently, an asymptotic expression for the correlation 
length has been derived for this case including the exact 
value of the prefactor 0]: 



e c / 2np, 

■ exp ( ^ 



8 2Trps 



1 - 



T 
Anps 



O 



Pi 



, (8) 



where e is the natural logarithm base. Here ^ is defined 
such that the dominant behavior of the pair correlation 
function at large distances is exp(— r/^) with power-law 
corrections. 

The quantum critical (QC) regime is the high tempera- 
ture scaling regime in the phase diagram of the u-model, 
derived by Chakravarty, Halperin and Nelson j^], where 
the dominant energy scale is set by the temperature: 



w - T 



(9) 



The crossover between RC and QC regimes occurs at 
T ~ const X ps ■ 

As the temperature increases even further, the assump- 
tions behind the mapping of the Heisenberg model to the 
(T-niodel eventually become invalid, which at tempera- 
tures T > c/a would destroy scaHng altogether. 

Chubukov and Sachdev recently argued that 
crossover from RC to QC regime occurs around T — 0.4 J 
for the S=l/2 model. According to their paper, in the 



region of T = 0.4J — 0.6 J, the system exhibits QC be- 
havior. The conclusions in ^ were based on the agree- 
ment of quantum critical scaling predictions for the bulk 
susceptibility, and the NMR spin-lattice relaxation 
rate, l/Ti, with the experimental measurements in 
La2Cu04 and with numerical calculations for the Heisen- 
berg model. 

However, the temperature dependence of the corre- 
lation length did not appear to support this picture of 
the QC behavior at intermediate temperatures ^(T) 
has been measured by means of inelastic neutron scat- 
tering in both La2Cu04 and Sr2Cu02Cl2 and cal- 
culated numerically for the 2D Heisenberg model using 
Monte-Carlo simulations ||l^ and series expansions ap- 
proaches |11|. All experimental measurements and nu- 



merical results agree well with each other, and the RC 
expression seems to describe well the combined experi- 
mental and numerical data for the correlation length not 
only for T < ||^ , but even for temperatures as high 
as T = J ^ ^Ps At T = J, the second term in 
the expansion in powers of T/ p^ is only about two times 
smaller than the first term, therefore Eq.(^ may remain 
accurate to all orders only if the coefficients in front of 
higher terms in the series in T / ps are anomalously small. 
Since these coefficients are presently unknown and may 
indeed turn out to be small, on the basis of the correla- 
tion length data alone one would be tempted to conclude 
that the system remains in the RC regime for tempera- 
tures up to T = J, which would also imply an absence of 
a QC regime for this model 

The main thesis of the present paper is that such a 
conclusion is not justified if one examines the detailed 
wave-vector and temperature dependences of 5'(q) and 
x(q). It has been shown in our earlier paper ||l^ that the 
scaling functions for the static susceptibility x(q)j ^ind 
equal time correlation function, S{c\), in the temperature 
range T — 0.6 J — J are that of the a-model in the QC 
rather than the RC regime. In that temperature range, 
we also obtained various universal dimensionless ratios in 
excellent agreement with their respective QC predictions, 
with the RC predictions being far outside the errorbars. 

However, the temperature dependence of the correla- 
tion length still presented a puzzle. In the quantum crit- 
ical region, the leading temperature dependence of the 
correlation length is of the form 



C = A, 



T' 



(10) 



where c is the spin-wave velocity and Xao is a univer- 
sal constant. With the known value of X^o |,|,|l| and 
the known zero temperature value of the spin- wave veloc- 
ity, the temperature dependence of the correlation length 
did not fit Eq.® in any extended temperature region, 
although the magnitude of ^ around T ~ J/2 is in agree- 
ment with Eq. (p^ . In Ref. ||l2|, we conjectured that 
there are nonuniversal lattice corrections at intermediate 
temperatures, which affect such quantities as the corre- 
lation length, but nevertheless the scaling functions and 



3 



the characteristic frequency of spin fluctuations are that 
of the QC regime. 

In the present paper, we are able to show that the 
disagreement of the correlation length with the quan- 
tum critical prediction is indeed a result of lattice cor- 
rections, more specifically corrections to the spin wave 
velocity, which is temperature dependent especially at 
higher temperatures. We further show that when Eq.(|l0|) 
is evaluated using the correct c(T) rather than its T = 
value, the agreement with the QC behavior is observed 
for temperatures above T ~ 0.55 J. 



III. UNIVERSAL SCALING FUNCTIONS 

The first stage of our analysis is to establish whether 
or not, and in what temperature range, universal q~ and 
cij— dependences of x(q, oj), derived for the quantum non- 
linear fj-model, apply to the S=l/2 and S—1 square lat- 
tice Heisenberg antiferromagnets. Here, it is convenient 
to introduce the spin wave mass defined as 
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m(T) 



c{T) 



(11) 



which, in what follows, replaces ^(T) as a scaling param- 
eter. We emphasize that both c{T) and £,{T) are mea- 
surable thermodynamic quantities. This replacement fa- 
cilitates scaling comparisons: we will show that when at 
temperatures T JS the spin wave velocity becomes 
temperature dependent due to lattice corrections to c, 
the behavior of m remains universal. 

For the u-model calculations, we employ expan- 
sion approach and utilize some of the results obtained 
earlier by Chubukov and Sachdev |^. We do not calcu- 
late the frequency dependence of the dynamical suscep- 
tibility explicitly. Instead, we base our comparisons on 
the calculated difference between the equal time correla- 
tors and static susceptibilities at all wavevectors, which 
is a sensitive probe of the frequency distribution of the 
spectral weight. Our calculations are discussed in detail 
in Appendix |^, where we derive the following equations 
which relate S{q) and x(q), computable via series expan- 
sions for the Heisenberg models, to the scaling parame- 
ters m and c: 



x(Q) 
x(q) 



5(q) 



1 



(-) 



Ql 



rx(q) 



where 

U2 ^ l + 0.1473/iV= 1.0491 
and the hmction 7 [a;] is defined according to: 



(12) 



(13) 



(14) 



x/2 



tanh(2:/2) 



for X > 1. 



(15) 



Here the expression for 1*2 is evaluated at TV = 3, which 
corresponds to the physical three-component spin. As 
discussed in Appendix these equations are not exact. 
In particular, there are nonzero higher order terms in 
|q— and multiplicative factors in front of c. However, 
such corrections are found to be less than 1 — 2%, which 
greatly simplifies further considerations. 

Eqs.( [l2| , [l3|) are in a form that greatly facilitates com- 
parison with the numerical data: the left hand parts 
of the above equations are computable numerically for 
the Heisenberg models, and their linear dependence on 
|q — Qp, expected from Eqs.( p^ , p^ ), can be verified by 
plotting on the appropriate scale. In the temperature 
range where the expected linear behavior is indeed ob- 
served, the scaling parameters m and c can be determined 
independently from the intercept and slope of the plotted 
data. 

Before we proceed with the comparisons, an important 
remark is in order. Since at high temperatures x(l) ^'^d 
S{q) do not exhibit strong g-dependence, their small de- 
viations from the values at Q = (tt, tt) are always linear 
in |q— Qp, as required by the lattice symmetry, so that at 
high enough temperatures q-dependences in Eqs.(^2|jl^) 
would convey no information about scaling. The high 
temperature crossover to non-universal weakly interact- 
ing spins is discussed more thoroughly in the next section. 

We start with Eq.(|l2|) and plot its left hand side, cal- 
culated using high temperature series expansions, versus 
|q — Qp for the S=l/2 and S=l Heisenberg models in 
Fig.|l|. Eq.(|T^) predicts that the plotted data should be 
linear. For wavevectors |q — Q| < 0.8, this is indeed 
the case for both S=l/2 and S=l models in a broad 
range of temperatures: T = 0.4 J — J for S=l/2, and 
T = 1.5J- 1.9J for S=l. 

We now turn to Eq.(^3|) and plot its left hand side 
versus |q— Qp in Fig||. Again, as expected from Eq.(p^, 
the plotted data is linear in the same wavevector and 
temperature range for both S=l/2 and S—1 models. 

Both Figsj^JI are generated with q — Q pointing along 
the diagonal of the Brillouin zone. The data with q — Q 
pointing along other directions is virtually identical in 
the same wavevector and temperature range, as expected 
in the isotropic a-model. Further, a straightforward esti- 
mate shows that corrections to linear spin wave spectrum 
due to higher gradients, not included in the a-model, 
are only about 5% for the wavevector range of interest, 
|q — Ql < 0.8. The studied temperature range is limited 
from below by our ability to evaluate accurately x(q) ^-nd 
5(q) from our high temperature series, presently known 
to /3" for the S=l/2 model, and to f3^° for the S=l 
model. 

From the slope and intercept of the linear fits to the 
data plotted in Fig.^, we are able to determine c and m 
separately. Since this data is used extensively in what 
follows, we carry out comparisons of our results for m 
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and c with direct calculations of the correlation length. 
The comparisons, which are described in Appendix 
confirm that our calculations are indeed accurate. 

We have thereby established that Eqs.( |l^ , |l3|) are valid 
in the range T = 0.4J - J {T ^ 1.5 J - 1.9J) for S=l/2 
(S=l) for |q — Ql < 0.8, which means that wavevector 
dependences of S{q) and x(q) in the Heisenberg models 
are described by the cr-model in the same wavevector and 
temperature range, which is shaded in Fig.|^. 



IV. SCALING REGIMES 

A. Quantum Nonlinear cr-Model 

We now proceed with the discussion of scaling regimes. 
The dynamical scaling form of the staggered correlator 
is 0: 



x.(g,c.) = x.(0)$«p,!^), 

V TO Uj/ 



(16) 



where the scaling function for the renormalized classical 
regime, differs from that in the quantum critical 

regime, ^Q*-'. In general, u> and m are not equivalent: 
the spin wave velocity, and therefore to, are defined at 
q ^ C^^, while u) is defined at g ^ S,^^- In the quantum 
nonlinear cr-model, c is expected to be temperature de- 
pendent only when the temperature is comparable to the 
ultraviolet cutoff A (except in case of relativistic invari- 
ance of the cutoff procedure, when c is always tempera- 
ture independent and coincides with its bare value). 

In the RC regime temperature dependences of the spin 
wave mass for = cxd and A^ = 3 are different: |,|: 



TO ~ exp 



T 



T forA^ = oo 
p. forA^^S ■ (1^) 



Furthermore, in this regime lo is not equivalent to to, 
except for N — oo: 



TO 



1 for A^ = oo 

(T/p,)i/2 forAr = 3 



(18) 



The exact asymptotic expression for m follows immedi- 
ately from Eq.(P): 



TO 



■ exp 



T 



T 



O 



J.2 



(19) 



The value of p^, which is the only dimensional parame- 
ter in Eq.(p^), is known for the Heisenberg models from 
zero temperature series expansions []l^ , quantum Monte 
Carlo simulations ||l^] , and l/S calculations . 

In the QC regime, temperature is the only energy scale 
in the problem, so that: 



TO ^ r. 



(20) 



which is valid for both A' = 3 and N = oo models. Fur- 
thermore, the ratio to/T is universal in the cr-model, and 
has been calculated via the 1 /N approach ||] : 



T 



= e 1 



0.2373 ^ 



0.9624 



0.2284 

A^ 



1.04, (21) 



where 9 = 2 log[(l + \/5)/2] . Given that already the first 
subleading 1/A^ correction is small (~ 8%) for A^ = 3, 
Eq.(pTl) is likely to be close to the true value for the 
0(3) cr-model, and is hereafter adopted as the QC pre- 
diction. It can be compared with calculations based on 
other techniques: 



T 



0.9, 

1.25 ± .25, 



2 + e expansions [3] 
Monte Carlo 111 



(22) 



The results are generally consistent with each other. Fur- 
ther support for the numerical value obtained in the 1/A^ 
calculation comes from the study of the two-plane model 
at its r = critical coupling by Sandvik and Scalapino 
1^, where temperature dependence of the correlation 
length was found to be in agreement with that derived 
from Eq.(|2^). 

Thus, both RC (^ and QC (|2l|) predictions for m{T) 
are available to carry out comparisons with Heisenberg 
models without any adjustable parameters. 



B. Heisenberg Models 

We start by comparing behavior of the spin wave mass 
for S=l/2, S=l, and the two plane models. The S—1 
model, due to its larger spin, is "more classical" than the 
S=l/2 model, hence it is less likely to have a range of 
quantum critical behavior. On the other hand, for the 
two plane model at its critical point, temperature depen- 
dences of xo and ^ consistent with the QC regime had 
been reported earlier by Sandvik and Scalapino ||^ . The 
comparison of different models proves quite informative. 

We now plot m{T) for these three models along with 
the universal RC ( [l9| ) and QC (21) scaling predictions 
in Fig.^ It is evident from Fig.^ that m{T) for the S=l 
model shows RC behavior at all temperatures, for the 
two plane model QC behavior at all temperatures, and 
for the S=l/2 single plane model a crossover between 
RC behavior below T ^ 0.45J and QC behavior above 
T = 0.65J. 

We now discuss the onset of nonuniversal behavior 
at high temperatures. The temperature scale at which 
nonuniversal behavior arises is set roughly by c/a. The 

and 



spin wave spectrum becomes nonlinear at q ^ ( 
the excitations in this wavevector range are suppressed 
only when T <C c/a. The spin stiffness, on the other 
hand, is less relevant to this crossover to nonuniversal be- 
havior at high temperatures. To a crude approximation 
we expect scaling to fail for S=l/2 and S=l models at 
T / (c/a) of the same order, given their similar spin wave 
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spectrum. In order to study failure of scaling, we there- 
fore plot the ratio m/T versus T/{co/a) (Fig.^. The 
use of Co instead of the temperature dependent c(T) is to 
avoid a nonlinear temperature scale. In this particular 
case, the difference between c and cq is less important 
given the ambiguous character of the upper boundary of 
scaling behavior. 

One can see from Fig.^ that deviations from the uni- 
versal behavior of m occurs at similar values of the ratio 
T/{co/a) ~ 0.6 for both S=l/2 and S=l models, above 
which a crossover to the T ^ J asymptotics m ~ 1/T oc- 
curs. The region where m becomes nonuniversal and no 
longer represents the spin wave mass is shaded in Fig||. 
Expressed in terms of J, the range where universal be- 
havior of m{T) holds is below T ~ J for the S=l/2 model 
and below T ^ 1.9J for the S=l model. 

Now, we turn to the analysis of the spin wave velocity 
calculated in Sect. Ill, and plot c determined from the 
data at different temperatures using Eq.|l| in Fig.| 

Dashed lines show the known T — value of the spin 
wave velocity calculated using 1/S expansions pq], 



= 2V2JS (1 + 0.1580/25' + 0.0216/(25)2) 

1.67 J forS=l/2 
3.07J for S=l ' 



(23) 



which for S=l/2 compares favorably with the T=0 series 
expansion calculations Q . 

For all temperatures studied the spin wave velocity 
turns out to be temperature dependent for both S=l/2 
and S=l models. Towards the lowest temperatures in 
the numerically accessible range, c approaches the ex- 
pected T=0 values given by Eq.(p3|). Since in the tr- 
model the spin wave velocity is temperature independent 
except near the ultraviolet cutoff, which in our case is 
set by the scale of order c/a ~ 2^/2JS, the observed 
temperature dependent spin wave velocity clearly results 
from the quantum nature of spin, when the temperature 
becomes comparable to JS. 

The quantum corrections to the spin wave velocity 
originating from higher order 1/5 terms in the spin wave 
theory, were calculated by Kaganov and Chubukov fl^ : 



c{T) = Co 



l + ^(- 
87r5 ' 



T 
SJ 



(24) 



where Ci^) is the Riemann zeta function. The predictions 
of Eq.(24), evaluated for S=l/2 and S=l, are shown as 
solid lines in Fig.^,b, respectively. The agreement with 
numerical calculations is quite good. Despite higher or- 
der terms in the expansion in powers of T/ JS, it appears 
that Eq.(24) describes the observed temperature depen- 
dence of c fairly well for both 5=1/2 and 5=1. 

At this stage an important question is, does the tem- 
perature dependence of c destroy scaling altogether? One 
might expect that once lattice corrections become im- 
portant for quantities like the spin-wave velocity, they 



might also lead to a breakdown of scaling altogether by 
introducing deviations from linearity in the spin wave 
spectrum and appearance of spectral weight in modes 
not contained in the quantum nonlinear cr-model. We 
argue that this is not the case here. Indeed, from the 
fact that Eqs.(|l^,^3|) remain valid when c(T) is already 
temperature dependent, we conclude that the spin wave 
spectrum still remains linear even though c{T) deviates 
from its T=0 value. By the same token, we conclude that 
there are no modes with significant spectral weight other 
than those described by the cr-model. 

Another universal quantity which has not been dis- 
cussed in our paper so far is the overall prefactor, xq, 
in front of the scaling expression for x(q, w) in Eq.(nq). 
This prefactor, which is renormalized by short wave- 
length quantum fluctuations, can be expressed in terms 
of iVo and ps for temperatures where corrections of or- 
der T/ JS can be ignored. Where c is T-dependent, this 
quantity becomes T-dependent as well. However, quan- 
tum corrections to the prefactor due to 1/5 terms have 
not yet been calculated analytically. 

Now the reason why the quantum critical scaling fails 
to describe the temperature dependence of the correla- 
tion length is clear. The correlation length derived from 
Eq. ([A^) is ^ = c/m. If one assumes that the spin wave 
velocity is temperature independent above T = 0.4 J, the 
quantum critical expression does not agree with the nu- 
merical data. Given the discussions above, we argue that 
one should consider not the temperature dependence of 
^, but that of m = c/^, because m is unaffected by the 
lattice corrections to c. Then, the crossover from RC 
behavior to QC behavior at temperatures around J/2 
becomes evident. 



V. CONCLUSION 

Magnetic properties of the doped superconducting 
cuprates may be closely related to the undoped parent 
compounds, such as La2Cu04 and Sr2Cu02Cl2, which 
are well described by the square lattice Heisenberg model. 
Motivated by this relationship, we have studied scaling 
and the role of lattice corrections in two dimensional 
square lattice Heisenberg antiferromagnets. 

In 1^, Chakravarty, Halperin, and Nelson mapped 2D 
coUinear quantum antiferromagnets onto the quantum 
nonlinear cr-model, and derived the corresponding phase 
diagram. When the T=0 ground state of the correspond- 
ing cr-model has Neel order, as is the case for the square 
lattice Heisenberg models, the low temperature scaling 
regime is classical (often called renormalized classical, or 
RC), where the characteristic energy scale ui ^T, while 
at higher temperatures a crossover to the quantum crit- 
ical (QC) regime may be observed, where to T. The 
S=l/2 square lattice Heisenberg model, which describes 
La2Cu04, was shown to exhibit RC behavior at low tem- 
peratures As it has a substantial ordered moment. 
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it was earlier believed not to have a region of quantum- 
critical behavior at intermediate temperatures. 

More recently, Chubukov and Sachdev ||] argued that 
the crossover from the RC to nonuniversal microscopic 
behavior occurs through the intermediate QC regime at 
0.4J < T < 0.6J, with universal behavior extending up 
to T = 0.6J and nonuniversal behavior above that tem- 
perature. The importance of this at first sight minor 
difference concerning a fairly narrow temperature range, 
arises from the fact that upon doping the QC region 
was argued to expand rather rapidly, covering a much 
broader, and well accessible temperature range 1). This 
scenario has important implications for our views of mag- 
netism in the superconducting cuprate materials ■ 

In our earlier paper we showed that the scaling 
functions for the equal time spin correlator, S'(q), and 
the static susceptibility, x{q), in the S=l/2 Heisenberg 
model in the temperature range T = 0.6 J — J agree with 
those calculated for the tr-model in the QC regime. The 
same turned out to be true for the ratio S (q) /Txiq) ■ Our 
comparisons did not involve any adjustable parameters. 

In the present paper, we undertook a detailed study of 
these issues using the high temperature series expansions 
approach. The simultaneous study of S=l/2 and S—1 
models allows us to contrast the two cases and highlight 
the parameter-free agreement of the intermediate tem- 
perature behavior for the 5" = 1/2 case with the universal 
quantum critical predictions of the non-linear cr-model. 
For both 5 = 1/2 and S = 1 models, we find that the 
q-dependences of S'(q) and x(q) separately, as well as of 
their ratio, agree very well up to T ~ 0.6co/a with the 
cr-model predictions calculated separately using 1/iV ex- 
pansions (for 1/iV calculations, we utilized many of the 
results obtained in [^). This allows us to calculate the 
temperature dependence of the spin wave mass, defined 
as TO = c/^, for which both RC and QC predictions are 
known without adjustable parameters. 

For S=l/2, we find that to(T) is in good agreement 
with the RC prediction for T < 0.45 J, while for 0.65 J < 
T < J the results are very close to the QC prediction 
and far from the RC prediction, exhibiting a crossover be- 
tween these two regimes for 0.45 J < T < 0.65 J. Further- 
more, we find that for T > 0.5J, the lattice corrections 
cause a temperature dependence of the spin-wave veloc- 
ity c and affect the overall prefactor in front of the scaling 
expression for xil,^)- The new results strongly support 
our earlier conjecture ||l2| that lattice corrections are such 
that full q- and w-dependences remain universal in the 
QC regime for temperatures up to T ~ 0.6co/a ~ J: 



c|q-Q| 
T 



(25) 



where <I'*^^(a;, y) is the scaling function of the QC regime, 
computable using, for instance, the approach of For 
the S=l model, where quantum fluctuations are substan- 
tially less, we find that RC behavior directly crosses over 
to the high temperature interacting local-spin behavior 



without an intervening QC regime. Finally, the observed 
lattice corrections to c for both S=l/2 and S=l roughly 
agree in magnitude and temperature dependence with 
earlier calculations of quantum T / JS effects by Kaganov 
and Chubukov [^9). 

In their recent publication, Grcven and coworkers 
pointed out the agreement of the measured and numeri- 
cally calculated correlation length ^(T) with RC predic- 
tion for all temperatures below T = J. In our work, we 
present an alternative theory which is not only consis- 
tent with the same data for £,{T), but explains the q- 
and a;-dependences of spin correlators as well. We show 
that detailed examination of the wavevector and temper- 
ature dependence of S{q) and x(g) supports the picture 
of renormalized classical behavior for temperatures up to 
T = 0.45 J, and a crossover to a quantum critical behav- 
ior around T = J/2. 

While the QC behavior for the spin wave mass occurs 
in our calculations in the temperature range T > 0.6J 
currently inaccessible to experiments in La2Cu04 and 
Sr2Cu02Cl2, many signatures of the QC behavior have 
been experimentally observed at lower temperatures. For 
instance, the nuclear relaxation rate, l/Ti, in La2Cu04 
saturates above T = 600K ~ 0.4J |7| at a value close 
in magnitude to the universal QC prediction This is 
consistent with the present study, as it has been shown 
|p|,pl| that for those quantities which in the RC regime 
depend on ^ through its logarithm only, such as the ra- 
tio T1T/T2G and the bulk susceptibihty, xo, the RC-QC 
crossover should be shifted substantially towards lower 
temperatures, in agreement with the measurements of 
Xo i and TiT/T2g § in La2Cu04. 

It has been proposed earlier that the temperature 
range of quantum critical behavior rapidly expands with 
doping, and that the scaling behavior described by the 
quantum nonlinear cr-model may be observed in doped 
antiferromagnets, and specifically in the high tempera- 
ture superconductors ||j2^,|l^. We hope that the study 
of the renormalized classical to quantum critical scaling 
crossover and of the origin of nonuniversal corrections to 
scaling in the insulator presented here, will be helpful in 
understanding magnetic behavior of the doped systems 
as well. 
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APPENDIX A: QUANTUM NONLINEAR 
a-MODEL 

The mapping [^ of the low energy spectrum of the 
Heisenberg model to the quantum nonlinear sigma model 
has been extensively discussed in the literature and text- 
books (e.g., see [Q). Here we briefly discuss some as- 
pects of this mapping, and of the expansion in this 
model, as necessary for our purposes. We refer the reader 
to Refs. [^,|| for a far more detailed treatment. 

The action of the quantum nonlinear cr-model can be 
written as: 



2 



l/T 



dT 



4 



[Vnf 



= 1, 



(Al) 



where the constraint = 1 captures the scattering ef- 
fects and gradient terms lead to the linear spin wave spec- 
trum. Here n is the antiferromagnetic order parameter, 
Pso is the bare spin stiffness and cq the bare spin wave ve- 
locity. The assumed ultraviolet cutoff is set by the lattice 
scale of the underlying Heisenberg model. 

This model has been extensively studied in recent years 
using rcnormalization group methods [pi, Monte Carlo 
simulations [^, and 1/A^ expansions The results 
obtained using different techniques agree with each other 
quite well. In our calculations, we rely mostly on the 
expansions approach. In the next Subsection, we discuss 
the leading order of this expansion, namely, the N = 
oo approximation, and then proceed with the subleading 
1/A^ corrections. 



1. N=oo Approximation 

The N — oo approximation is based on the replace- 
ment of the physical 0(3) quantum nonlinear cr-model, 
where the three-component order parameter n corre- 
sponds to three-component physical spin, by the 0(oo) 
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cr-model where the number of components of n is infinite. 
A consequence of such a replacement is the absence of 
damping. This makes the N — oo approximation unsuit- 
able for such problems as the nuclear spin-lattice relax- 
ation or the neutron scattering, where correct description 
of damping is essential. Nevertheless, for our particular 
purposes of extracting scaling parameters from the data, 
already the N = oo approximation is relatively accurate. 
It is further improved by taking into account the first 
subleading corrections, which are discussed in the next 
Subsection. 

The unit length constraint can be enforced using La- 
grange multiplier. A: 



l/T 



dT 



— {drfif + {Vnf + i\ {n^ - 1) 



(A2) 



The action is now quadratic in n which can therefore 
be integrated out, leaving a functional integral for the 
field A. In case of iV = oo, the A field does not fluctuate 
around the saddle point value i(A), which in what follows 
is denoted as for the reason which becomes clear mo- 



mentarily. By performing Fourier transform of Eq.(A2) 
one obtains the N = oo solution for the Matsubara fre- 
quencies H: 



A 



c 

TO ' 



(A3) 



where c is the spin wave velocity and it is now evident 
that TO = c/^ has the meaning of the spin wave mass, i.e. 
the gap in the spin wave spectrum. 

The Matsubara correlator Xs(g, iaj^), where Un — 
27rnT, is related to the physical real frequency response 
function by analytical continuation: uj iun- It is easy 
to see that Eq.(A3) corresponds to a spectrum of un- 
damped magnons with a gap: 



A 



TO^ -I- c'q 



,2„2 



A 
2^ 



1 



1 



where the spin wave dispersion is: 



(A4) 



(A5) 



The imaginary part of which describes dissipa- 

tion, is: 

x'^{q,Lj) = — [S{u;~eq)-S{cj + eg)] at N^oo. 

(A6) 

Of course, for the physical 0(3) QNLcr model the damp- 
ing is finite and the peaks in x" around oj = ±£q acquire 



finite width. This effect is captured only in the terms 
beyond the N = oo approximation. 

Our considerations for the S=l/2 and S—1 Heisenberg 
models are based on the high temperature series data 
for the q-dependent static susceptibility of the staggered 
order parameter (n) field , Xsil), and equal time struc- 
ture factor of the same field, S{q). Using Eq.(A3), one 
obtains: 



x{q) = Xs{q,i^n = 0) = ^, 



+ 00 



A 



tanh ^ 



n— — oo 



2T 



where 



A^XsiO) X TO^ 



(to^ + c^q'^) 



It follows directly from Eq.(lA7|) that: 
Xs{0) , 



Xs{q) 



q for TV = oo. 



(A7) 



(A8) 



(A9) 



(AlO) 



Further, the ratio of Eq.(A7) and Eq.(A8) depends only 
on Eq/T: 



S{q) 



^tanh-i 



2T 



TXs{q) 2T 
We now introduce a function ^[x\ defined as: 
x/2 



for a; > 1, 



(All) 



(A12) 



tanh(a::/2 

which has the following asymptotic behaviors: ^[x\ ~ 
2^3(2: - 1) for < a; - 1 < 1 and ^[x] ~ 2x for a; > 1. 
Then 



S{q) 



TXsiq) 



J, J ■ yj;) 9' for TV = 00. (A13) 



2. 1/iV Expansion 



In this Subsection, we discuss corrections to the N 



00 equations ( AlC , Al^ ) for the physical case of = 3. 
We utilize the results of the 1/iV expansion calcula- 
tions by Chubukov, Sachdev, and Ye 1^ and find that 
Eqs.( |Al^ , |ATo| ) remai n quit e accurate even for iV = 3. 

We turn first to Eq.( A13 ). While the full q-dependence 
of 1/iV corrections have not been calculated so far, the 
results for small wavevectors as presented in |^ turn 
out to be quite sufficient for our purposes. In the 
asymptotic high temperature (quantum critical) regime, 
to/T is temperature independent, while XsiQ)/Xsiq) and 
S{q)/Txs{q) depend only on the product q£, ~ cq/T. 
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With straightforward appHcation of the results of |^ , we 
obtain for small q: 



been carried o ut. Instead, in the next Section we verify 



XsjO) 
Xsiq) 



= 1 



vi— q 

m/ 



S{q) 



TXsiq) 



(A14) 



(A15) 



where 



vi = \~ 0.0042/iV-^ 0.9986, 
W2 = 1 + 0.1473/iV^ 1.0491, 
V2 = 1 + 0.0382/iV^ 1.0127 



(A16) 



(1/iV terms are evaluated for = 3). 

The 1/N corrections to the coefhcients vi, V2, which 



enter in parenthesis in front of q^ in Eqs.( A14 A15 ), are 
essentially negligible. Higher terms in the expansion in 
powers of q"^ have not been calculated so far. However, 
once the leading terms have very small corrections, cor- 
rections to higher order terms in the expansion in 
should be of the same order of magni tude. Th e reason 
for this conclusion is as follows. Eqs.(Al4, A15) can be 



considered as relationshps between the position of the 
closest to the origin singuliarity for imaginary wavevec- 
tors, gsing = 1/C) and the expansion in powers of q^ near 
q = 0. Were the higher order terms large, the position 
of the singuliarity expressed in terms of the coefhcient 
in front of g^, namely, vi or V2, would change, which in 
its turn would lead to large corrections to vi,V2- Given 
the fact that such corrections are in fact small, and ig- 
noring the possibility of accidental cancellations for both 
Eqs.( A14, Al£ ), we arrive to the conclusion t hat highe r 
order g^" terms are small and therefore Eqs.( A14, Al£ ) 
are applicable for ga <C 1 , not only for the much narrower 
wavector range <C 1. 

It is quite remarkable that the 1/N correction to the 
coefficient U2 is small even when evaluated at iV = 3. 
Indeed, the ratio S{Q)/Tx{Q) reflects the distribution 
of the spectral weight over frequencies. At the com- 
mensurate wavevector, the character of this distribution 
changes from a delta- function peak at finite w = m to a 
smooth peak centered at w = as one goes from N = oo 
to N = 3, which in general could have led to large 1/iV 
corrections to U2. The fact that the rigorously calcu- 



lated correction (A16) is in fact small shows that already 
the N — oo approximation accurately captures the re- 
lationship between the spin wave mass and the relative 
correlation strength at different Matsubara frequencies. 

Further, since the above corrections are essentially neg- 
ligible in the QC regime where they have been calcu- 
lated, it iseems reasonable to conjecture that they may 
not increase drastically as m/T decreases, in which case 
they should apply to the RC-QC crossover regime as 
well. The rigorous calculations to check this conjec- 
ture are prohibitively complicated and therefore have not 



Eqs.( A14 , A15 ) are in agreement with the numerical data 
for Hcisenberg models for all temperatures studied, and 
thereby establish that the above conjecture is indeed cor- 
rect, which greatly simplifies futhcr considerations. 

To summarize, in this Appendix we derived surpris- 
ingly simple Eqs.( A14, |A15| , Al^ ). Negligible (for vi, V2) 
or small (~ 5% for U2) 1/N corrections to the N — 00 
result indicate a well behaved large N expansion, and 
we expect the above equations to be accurate for the 
physical 0(3) model. The correspondence between the 
Heisenberg and sigma-models is such that Xs{q,^) for 
the cr-model is x(|q — QL^^) for the Heisenberg model, 
wher e Q = (tt , tt). Taking this into account, we rewrite 



Eqs.( A14, A15) in the form appropriate for the Hcisen- 
berg models, Eqs.( p"2|jr^ ), and use them in Section HI for 
comparisons with the numerical data. 



APPENDIX B: CORRELATION LENGTH 

In this Appendix, we verify our resluts for m and c 
by comparing the ratio c/m with the independently cal- 
culated correlation length. We calculate the correlation 
length, defined as the inverse rate of exponential decay of 
spin correlator at large distances, directly using our new 
"imaginary wavevector" method |Q. The idea of this 
calculation is based on an observation that when the pair 
correlation function decays exponentially up to power- 
law corrections, which is the case for two-dimensional 
Heisenberg models, the Fourier transform of spin corre- 
lator around the ordering wavevector, x(Q + 11^(7), where 
rig is a unit vector, should have a singuliarity on the 
imaginary q axis at Qsing — i/S, when analytically contin- 
ued to complex q. The position of this singuliarity can 
be determined through high temperature series expan- 
sions using the methods described in |2^. In Fig.|^, our 
results for ^ for the S=l/2 Heisenberg model are plotted 
along with the earlier quantum Monte Carlo calculations 
by Makivic and Ding the agreement is excellent at 
all temperatures. 

Now we can verify the equality ^ = c/m, which holds 
by the definition of m. In Fig.|[ we plot our directly cal- 
culated ^ by a solid line, the ratio c/m determined solely 
from Eq.(|2|) as hollow circles, and the ratio of c and m 
with both determined seperately from Eq.(13) as solid 
circles. The agreement between all three sets of data is 
excellent for both S=l and S=l/2. This agreement not 
only indicates that our calculations of m and c are ac- 
curate, but also serves as another verification that the 
cr-model corectly describes wavevector and frequency de- 
pendences of the dynamical susceptibility of the Heisen- 
berg models in the studied temperature range. 



FIG. 1. x(Q)/x(ci) is plotted as a function of |q— Q|^ for 
the S=l/2 (a) and S—l (b) Heisenberg models for selected 
temperatures. Linear dependences, in agreement with those 
derived for the cr-model Eq.(^2|), are observed. 
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FIG. 2. 7^ [5'(Q)/rx(q)] is plotted as a function of 
|q — Ql^ for the S=l/2 (a,aa) and S=l (b) Heisenberg mod- 
els for selected temperatures. The dependence is linear, in 
agreement with those derived for the a-model Eq.(p^. For 
the S=l/2 model the intercept at |q — Q| = is nearly tem- 
perature independent above T = 0.6J, and is in agreement 
with the known universal QC prediction ~ 1.19, shown as an 
arrow. For temperatures lower than T — 0.5 J for S=l/2, and 
for all temperatures for S=l, the intercept is smaller than the 
QC prediction and decreases as the temperature decreases, 
indicating crossover to RC behavior. 

FIG. 3. A drawing of the Brillouin zone for the square 
lattice; shaded area shows the wavevector range where spin 
correlations of the Heisenberg models are described by the 
CT-model. 



FIG. 7. The correlation length determined using our new 
method of analytical continuation to imaginary wavevectors 
(solid line) plotted along with earlier Monte Carlo results or 
Makivic and Ding ||l^ (circles). 

FIG. 8. The correlation length for the S=l Heisenberg 
model calculated directly using the analytical continuation to 
imaginary wavevectors (solid line), and independently u sing 
the mapping to the quantum nonlinear cr-model, from Eq.(|l2) 
(triangles) and Eq.(|l3|) (circles). All three sets of data agree 
with each other for both S=l/2 and S=l. 



FIG. 4. The ratio m/T is plotted for (a) the S=l/2, (b) 
S=l and (c) the two-plane Heisenberg models versus T/J 
along with the known RC (Eq.(p^), dot-dashed line) and QC 
(Eq.(pl|), dashed line) universal predictions. The compar- 
isons do not contain any adjustable parameters, because the 
spin stiffness entering the RC expression is known from earlier 
studies [ p"^Jl^ . For the "more classical" S=l model renormal- 
ized classical behavior, and for the "critical" two plane model 
quantum critical behavior, are observed at all temperatures. 
In contrast, for the S=l/2 model a crossover from RC to QC 
behavior occurs around T ~ 0.55 J. 



FIG. 5. m/T is plotted versus the ratio T/{co/a) for the 
S=l/2 (circles) and S=l (diamonds) single-plane models, as 
well as the model of two S=l/2 planes (triangles) coupled 
such that the system is at the critical point (as determined 
by Sandvik and Scalapino) and therefore should exhibit QC 
behavior at all temperatures. The difference between mod- 
els in the temperature range where scaling applies is evident: 
the S—1 model is always renormalized classical, the two plane 
model is always quantum critical, and the S=l/2 model ex- 
hibits a crossover between the two around T — 0.5J. More- 
over, for S=l/2 and S=l models m/T deviates from the re- 
spective scaling temperature dependences at roughly the same 
value of T/{co/a) > 0.6 (shown as shaded area), as expected 
given the similar character of deviations from linearity in their 
spin wave spectrum (see text). For co, we use the spin wave 
theory results: co = 1.67J for S=l/2, cq = 3.07J for S=l 
[^, and Co = 1.9J for the two-plane model [ |l7t . 



FIG. 6. The spin wave velocity plotted as a function of 
temperature for the S=l/2 (a) and S=l (b) Heisenberg mod- 
els. Dashed line shows the known T=0 values, co [ p^ , p^ , 
and solid line is the result of 1/S calculation by Kaganov 
and Chubukov [^, Eq.^. The temperature dependence for 
T > JS arises from nonuniversal lattice corrections. 
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I. INTRODUCTION 

The problem of antiferromagnetism in two dimen- 
sions has recently attracted much attention, in part 
due to its relevance to high temperature supercon- 
ductivity. Since the high temperature superconduct- 
ing materials are derivatives, obtained by doping, of 
parent insulating layered S=l/2 antiferromagnets, the- 
ory of the insulating phase is important for under- 
standing the unusual magnetic and transport proper- 
ties in the normal state of doped materials. Our re- 
sults are directly relevant to the following two materi- 
als: La2Cu04, which is the insulating parent compound 
of the La2-xSrxCu04 and La2-xBaxCu04 superconduc- 
tor families, and Sr2Cu02Cl2, which is unrelated to any 
superconducting family, but has well studied magnetic 
properties, similar to those of La2Cu04. Both materi- 
als experience antiferromagnetic transitions to the Neel 
state with nonzero staggered magnetization, La2Cu04 at 
Tn ~ 300 A' and Sr2Cu02Cl2 at Tn ~ 256 A'. According 
to the Hohenberg-Mermin- Wagner theorem, in genuine 
two dimensional antiferromagnets the long range order at 
finite temperatures is prevented by thermal fiuctuations. 
In real materials, small interlayer exchange interaction al- 
lows the transition to occur at a finite temperature Tjy. 
In what follows we discuss only the properties well above 
Tjv, where the magnetic subsystem can be considered 
two-dimensional. 

It is well established [1] that the magnetic properties of 
La2Cu04 and Sr2Cu02Cl2 are described by the 5=1/2 
Heisenberg model on the two-dimensional square lattice. 
The Heisenberg model is defined by the following Hamil- 
tonian: 



F = J^S,S,, (1) 

where Si are local spin operators and the notation ^^(jj) 
means that the sum is taken over all pairs of nearest 
neighbors. In what follows we set ks = Ti = a = 1, where 
a is the lattice constant, although in certain expressions 
we retain a for clarity. The dynamical spin susceptibility 
x(q,w) = x' + *x") is defined such that the magnetiza- 
tion induced by a wavevector and frequency dependent 
magnetic field is: 

M(q,c^) = x(q,w)Hext(q,c^). (2) 

We calculate high temperature series expansions in pow- 
ers of /3 = J/T for two quantities: the static q-dependent 
susceptibility, defined as a response to a static, but gen- 
erally nonuniform, magnetic field, 

x(q) = x(q,^ = 0), (3) 

and the Fourier transform of the equal time correlation 
function of two spin operators, 

5(q) = E(S-(o)S-W)e^p('q^)' (4) 

r 

where angular brackets represent thermal averages. 
These quantities can be expressed in terms of x"(q,'^), 
which plays the role of the spectral weight for spin fiuc- 
tuations, according to: 

x(q) = x(q, «t^n = 0) = - / duj (5) 

where w„ = 2'iTnT is the imaginary bosonic Matsubara 
frequency. 

Total spin conservation in the model leads to the rela- 
tion 

Xo = ffVljX(q = 0) = ffVl,T-i5(q = 0), (6) 

where xo is the bulk magnetic susceptibility, g the elec- 
tron gyromagnetic ratio, and /j^b the Bohr magneton. 

We study both the case of S = 1/2, which corresponds 
to La2Cu04 and Sr2Cu02Cl2, and the case of S = 1. 
We also study a model of two S=l/2 coupled Heisenberg 
planes at its critical point, J2/J1 — 2.51 ± 0.02 [2]. For 
both x(q) and S(q), we are able to generate series for 
an arbitrary wavevector q complete to order {3^^ for the 
5=1/2 model, and to l3^° for the S = I model. For 
the two-plane model we set J2/J1 = 2.5, and generate 
series for the structure factor and the susceptibility at the 
antiferromagnetic ordering vector to order l3^° . We use 
Fade approximation techniques for series analysis. We 
emphasize that the series converges extremely well in the 
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temperature region presented in the paper. In fact, with 
the available long series, most of the Pade approximants 
are within a per cent, or so, from each other, down to 
T ^ 0.4J for S = 1/2 and T ^ 1.5J for 5=1. This 
allows us to put quite small errorbars on our results for all 
but the lowest temperatures. At the lowest temperatures, 
the approximants begin to significantly differ from each 
other and hence we do not go to still lower temperatures. 

II. OVERVIEW OF RENORMALIZED 
CLASSICAL VERSUS QUANTUM CRITICAL 
SCALING BEHAVIOR 

Scaling behavior of the 2D Heisenberg models follows 
from their mapping [3] to the quantum nonlinear sigma 
model ((T-model hereafter). When the long range order 
is present at T = 0, fiuctuations at low enough tempera- 
tures are purely classical [3], in other words the character- 
istic frequency of spin fiuctuations Co <^T . This happens 
because for T Ps (Ps is the T=0 spin stiffness), 

~Ti/2exp(-27rp,/T) (7) 

decreases faster than the temperature as the tempera- 
ture decreases. Accordingly, this regime is called the 
renormahzed classical (RC) regime, where "renormal- 
ized" means that the T = parameters, such as the 
spin stiffness, ps, and spin wave velocity, c, are renormal- 
ized compared to their mean-field values as a result of 
quantum fiuctuations at shorter wavelengths. 

Recently, an asymptotic expression for the correlation 
length has been derived for this case including the exact 
value of the prefactor [4] : 




where e is the natural logarithm base. Here ^ is defined 
such that the dominant behavior of the pair correlation 
function at large distances is exp(— r/^) with power-law 
corrections. 

The quantum critical (QC) regime is the high tempera- 
ture scaling regime in the phase diagram of the cr-model, 
derived by Chakravarty, Halperin and Nelson [3], where 
the dominant energy scale is set by the temperature: 

w ~ T (9) 

The crossover between RC and QC regimes occurs at 
T ~ const X Ps ■ 

As the temperature increases even further, the assump- 
tions behind the mapping of the Heisenberg model to the 
(T-model eventually become invalid, which at tempera- 
tures T > c/a would destroy scaling altogether. 

Chubukov and Sachdev [5] recently argued that 
crossover from RC to QC regime occurs around T = 0.4 J 
for the S=l/2 model. According to their paper, in the 



region of T = 0.4J — 0.6J, the system exhibits QC be- 
havior. The conclusions in [5] were based on the agree- 
ment of quantum critical scaling predictions for the bulk 
susceptibility, xo, and the NMR spin-lattice relaxation 
rate, 1/Ti, with the experimental measurements [6,7] in 
La2Cu04 and with numerical calculations for the Heisen- 
berg model. 

However, the temperature dependence of the corre- 
lation length did not appear to support this picture of 
the QC behavior at intermediate temperatures [8]. C(T) 
has been measured by means of inelastic neutron scat- 
tering in both La2Cu04 [9] and Sr2Cu02Cl2 [8], and cal- 
culated numerically for the 2D Heisenberg model using 
Monte-Carlo simulations [10] and series expansions ap- 
proaches [11]. All experimental measurements and nu- 
merical results agree well with each other, and the RC 
expression seems to describe well the combined experi- 
mental and numerical data for the correlation length not 
only for T pg [3], but even for temperatures as high 
as T = J ~ 5/9s [8]. At T = J, the second term in 
the expansion in powers of T/ps is only about two times 
smaller than the first term, therefore Eq.(8) may remain 
accurate to all orders only if the coefficients in front of 
higher terms in the series in T/ ps are anomalously small. 
Since these coefficients are presently unknown and may 
indeed turn out to be small, on the basis of the correla- 
tion length data alone one would be tempted to conclude 
that the system remains in the RC regime for tempera- 
tures up to T = J, which would also imply an absence of 
a QC regime for this model [8]. 

The main thesis of the present paper is that such a 
conclusion is not justified if one examines the detailed 
wave-vector and temperature dependences of S(q) and 
x(q). It has been shown in our earlier paper [12] that the 
scaling functions for the static susceptibility x(ct)) and 
equal time correlation function, S(q), in the temperature 
range T = 0.6J — J are that of the cr-model in the QC 
rather than the RC regime. In that temperature range, 
we also obtained various universal dimensionless ratios in 
excellent agreement with their respective QC predictions, 
with the RC predictions being far outside the errorbars. 

However, the temperature dependence of the correla- 
tion length still presented a puzzle. In the quantum crit- 
ical region, the leading temperature dependence of the 
correlation length is of the form 

i = X-'^, (10) 

where c is the spin-wave velocity and X^o is a univer- 
sal constant. With the known value of X^o [3,5,13] and 
the known zero temperature value of the spin- wave veloc- 
ity, the temperature dependence of the correlation length 
did not fit Eq.(lO) in any extended temperature region, 
although the magnitude of ^ around T ~ J/2 is in agree- 
ment with Eq.(lO). In Ref. [12], we conjectured that 
there are nonuniversal lattice corrections at intermediate 
temperatures, which affect such quantities as the corre- 
lation length, but nevertheless the scaling functions and 
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the characteristic frequency of spin fluctuations are that 
of the QC regime. 

In the present paper, we are able to show that the 
disagreement of the correlation length with the quan- 
tum critical prediction is indeed a result of lattice cor- 
rections, more specifically corrections to the spin wave 
velocity, which is temperature dependent especially at 
higher temperatures. We further show that when Eq.(lO) 
is evaluated using the correct c(T) rather than its T = 
value, the agreement with the QC behavior is observed 
for temperatures above T ~ 0.55J. 



III. UNIVERSAL SCALING FUNCTIONS 

The first stage of our analysis is to establish whether 
or not, and in what temperature range, universal q— and 
w— dependences of x(q,i^), derived for the quantum non- 
linear (T-model, apply to the S=l/2 and S=l square lat- 
tice Heisenberg antiferromagnets. Here, it is convenient 
to introduce the spin wave mass defined as 
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m(T) 



c(T) 
e(T)' 



(11) 



which, in what follows, replaces C(T) as a scaling param- 
eter. We emphasize that both c(T) and C(T) are mea- 
surable thermodynamic quantities. This replacement fa- 
cilitates scaling comparisons: we will show that when at 
temperatures T ~ JS the spin wave velocity becomes 
temperature dependent due to lattice corrections to c, 
the behavior of m remains universal. 

For the cr-model calculations, we employ l/N expan- 
sion approach and utilize some of the results obtained 
earlier by Chubukov and Sachdev [5]. We do not calcu- 
late the frequency dependence of the dynamical suscep- 
tibility explicitly. Instead, we base our comparisons on 
the calculated difference between the equal time correla- 
tors and static susceptibilities at all wavevectors, which 
is a sensitive probe of the frequency distribution of the 
spectral weight. Our calculations are discussed in detail 
in Appendix A, where we derive the following equations 
which relate S(q) and x(ct)) computable via series expan- 
sions for the Heisenberg models, to the scaling parame- 
ters m and c: 
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x(Q) 
x(q) 



1 



TJ 



|q-Ql 



^)'|q-QP 



where 

U2 ~ l-F 0.1473/#= 1.0491 
and the function ')[x\ is defined according to: 



(12) 



(13) 



(14) 



x/2 



tanh(a;/2) 



for X > 1. 



(15) 



Here the expression for U2 is evaluated at # = 3, which 
corresponds to the physical three-component spin. As 
discussed in Appendix A, these equations are not exact. 
In particular, there are nonzero higher order terms in 
|q— Qp and multiplicative factors in front of c. However, 
such corrections are found to be less than 1 — 2%, which 
greatly simplifies further considerations. 

Eqs.(12,13) are in a form that greatly facilitates com- 
parison with the numerical data: the left hand parts 
of the above equations are computable numerically for 
the Heisenberg models, and their linear dependence on 
|q — Qp, expected from Eqs.(12,13), can be verified by 
plotting on the appropriate scale. In the temperature 
range where the expected linear behavior is indeed ob- 
served, the scaling parameters m and c can be determined 
independently from the intercept and slope of the plotted 
data. 

Before we proceed with the comparisons, an important 
remark is in order. Since at high temperatures x(q) ^nd 
^(q) do not exhibit strong (/-dependence, their small de- 
viations from the values at Q = (tt, tt) are always linear 
in |q— Qp, as required by the lattice symmetry, so that at 
high enough temperatures (/-dependences in Eqs.(12,13) 
would convey no information about scaling. The high 
temperature crossover to non-universal weakly interact- 
ing spins is discussed more thoroughly in the next section. 

We start with Eq.(12) and plot its left hand side, cal- 
culated using high temperature series expansions, versus 
|q — Qp for the S=l/2 and S=l Heisenberg models in 
Fig.l. Eq.(12) predicts that the plotted data should be 
linear. For wavevectors |q — Q| < 0.8, this is indeed 
the case for both S=l/2 and S=l models in a broad 
range of temperatures: T = 0.4 J — J for S=l/2, and 
T= 1.5J-1.9Jfor S=l. 

We now turn to Eq.(13) and plot its left hand side 
versus |q— Qp in Fig. 2. Again, as expected from Eq. (13), 
the plotted data is linear in the same wavevector and 
temperature range for both S=l/2 and S=l models. 

Both Figs. 1,2 are generated with q — Q pointing along 
the diagonal of the Brillouin zone. The data with q — Q 
pointing along other directions is virtually identical in 
the same wavevector and temperature range, as expected 
in the isotropic (T-model. Further, a straightforward esti- 
mate shows that corrections to linear spin wave spectrum 
due to higher gradients, not included in the (T-model, 
are only about 5% for the wavevector range of interest, 
|q — Ql < 0.8. The studied temperature range is limited 
from below by our ability to evaluate accurately x(q) ^nd 
^(q) from our high temperature series, presently known 
to for the S=l/2 model, and to (3^° for the S=l 
model. 

From the slope and intercept of the linear fits to the 
data plotted in Fig. 2, we are able to determine c and m 
separately. Since this data is used extensively in what 
follows, we carry out comparisons of our results for m 
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FIG. 1. x(Q)/x(q) is plotted as a function of |q — Q|^ for 
the S=l/2 (a) and S=l (b) Heisenberg models for selected 
temperatures. Linear dependences, in agreement with those 
derived for the cr-model Eq.(12), are observed. 

and c with direct calculations of the correlation length. 
The comparisons, which are described in Appendix B, 
confirm that our calculations are indeed accurate. 

We have thereby established that Eqs.(12,13) are valid 
in the range T = 0.4J - J (T = 1.5J - 1.9J) for S=l/2 
(S=l) for |q — Ql < 0.8, which means that wavevector 
dependences of S(q) and x(ct) in the Heisenberg models 
are described by the cr-model in the same wavevector and 
temperature range, which is shaded in Fig. 3. 



IV. SCALING REGIMES 

A. Quantum Nonlinear cr-Model 

We now proceed with the discussion of scaling regimes. 
The dynamical scaling form of the staggered correlator 
IS [3]: 



(0 (21 ^ 



(16) 



f [5(Q)/Tx(q)] = 

(«W^r + (c/T)^|q-QP 
Vf^oT 




0.2 0.4 0.6 O.S 
,T=0.4 

(aa) 




0.4 

iq-Qr 

FIG. 2. 7^ [5(Q)/Tx(q)] is plotted as a function of |q-Q|^ 
for the S=l/2 (a,aa) and S=l (b) Heisenberg models for se- 
lected temperatures. The dependence is linear, in agreement 
with those derived for the cr-model Eq.(13). For the S=l/2 
model the intercept at |q — Q| = is nearly temperature 
independent above T = 0.6J, and is in agreement with the 
known universal QC prediction ~ 1.19, shown as an arrow. 
For temperatures lower than T = 0.5 J for S=l/2, and for 
all temperatures for S=l, the intercept is smaller than the 
QC prediction and decreases as the temperature decreases, 
indicating crossover to RC behavior. 

where the scaling function for the renormalized classical 
regime, differs from that in the quantum critical 

regime, In general, ui and m are not equivalent: 

the spin wave velocity, and therefore m, are defined at 
while u) is defined at q ^ In the quantum 
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FIG. 3. A drawing of the Brillouin zone for the square 
lattice; shaded area shows the wavevector range where spin 
correlations of the Heisenberg models are described by the 
(7-model. 

nonlinear cr-model, c is expected to be temperature de- 
pendent only when the temperature is comparable to the 
ultraviolet cutoff A (except in case of relativistic invari- 
ance of the cutoff procedure, when c is always tempera- 
ture independent and coincides with its bare value). 

In the RC regime temperature dependences of the spin 
wave mass for N = co and N = 3 are different: [3,5]: 



exp 



T 



T for N = CO 
p, for # = 3 



(17) 



Furthermore, in this regime to is not equivalent to m, 
except for N = oo: 



U) 

m 



1 



iT/Ps) 



1/2 



for N 
for N 



oo 
3 



(18) 



The exact asymptotic expression for m follows immedi- 
ately from Eq.(8): 



■ exp 



-27rps 
T 



1 



T 



O 



(19) 



The value of pg , which is the only dimensional parame- 
ter in Eq.(19), is known for the Heisenberg models from 
zero temperature series expansions [14], quantum Monte 
Carlo simulations [15], and 1/ S calculations [16]. 

In the QC regime, temperature is the only energy scale 
in the problem, so that: 



(20) 



which is valid for both # = 3 and N = oo models. Fur- 
thermore, the ratio m/T is universal in the cr-model, and 
has been calculated via the l/N approach [5]: 



m 
T 



6(1 



0.2373 ^ 

N ' 



0.9624- 



0.2284 

N 



1.04, (21) 



where = 21og[(l-|-V5)/2]. Given that already the first 
subleading l/N correction is small (~ 8%) for # = 3, 
Eq.(21) is likely to be close to the true value for the 
0(3) (T-model, and is hereafter adopted as the QC pre- 
diction. It can be compared with calculations based on 
other techniques: 



T 



0.9, 

1.25 ± .25, 



2 -\- e expansions [3] 
Monte Carlo [13] 



(22) 



The results are generally consistent with each other. Fur- 
ther support for the numerical value obtained in the l/N 
calculation comes from the study of the two-plane model 
at its T = critical coupling by Sandvik and Scalapino 
[2], where temperature dependence of the correlation 
length was found to be in agreement with that derived 
from Eq.(21). 

Thus, both RC (19) and QC (21) predictions for m(T) 
are available to carry out comparisons with Heisenberg 
models without any adjustable parameters. 



B. Heisenberg Models 

We start by comparing behavior of the spin wave mass 
for S=l/2, S=l, and the two plane models. The S=l 
model, due to its larger spin, is "more classical" than the 
S=l/2 model, hence it is less likely to have a range of 
quantum critical behavior. On the other hand, for the 
two plane model at its critical point, temperature depen- 
dences of xo and ^ consistent with the QC regime had 
been reported earlier by Sandvik and Scalapino [2]. The 
comparison of different models proves quite informative. 

We now plot m(T) for these three models along with 
the universal RC (19) and QC (21) scaling predictions 
in Fig. 4. It is evident from Fig. 4 that m(T) for the S=l 
model shows RC behavior at all temperatures, for the 
two plane model QC behavior at all temperatures, and 
for the S=l/2 single plane model a crossover between 
RC behavior below T ~ 0.45J and QC behavior above 
T= 0.65J. 

We now discuss the onset of nonuniversal behavior 
at high temperatures. The temperature scale at which 
nonuniversal behavior arises is set roughly by c/a. The 
spin wave spectrum becomes nonlinear at q ~ a~^, and 
the excitations in this wavevector range are suppressed 
only when T ^ c/a. The spin stiffness, on the other 
hand, is less relevant to this crossover to nonuniversal be- 
havior at high temperatures. To a crude approximation 
we expect scaling to fail for S=l/2 and S=l models at 
T/(c/a) of the same order, given their similar spin wave 
spectrum. In order to study failure of scaling, we there- 
fore plot the ratio m/T versus T/(co/a) (Fig. 5). The 
use of Co instead of the temperature dependent c(T) is to 
avoid a nonlinear temperature scale. In this particular 
case, the difference between c and cq is less important 
given the ambiguous character of the upper boundary of 
scaling behavior. 
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FIG. 4. Tlie ratio m/T is plotted for (a) the S=l/2, (b) 
S=l and (c) the two-plane Heisenberg models versus T/J 
along with the known RC (Eq.(19), dot-dashed line) and QC 
(Eq.(21), dashed line) universal predictions. The compar- 
isons do not contain any adjustable parameters, because the 
spin stiffness entering the RC expression is known from earlier 
studies [14,16]. For the "more classical" S=l model renormal- 
ized classical behavior, and for the "critical" two plane model 
quantum critical behavior, are observed at all temperatures. 
In contrast, for the S=l/2 model a crossover from RC to QC 
behavior occurs around T ~ 0.55 J. 

One can see from Fig. 5 that deviations from the uni- 
versal behavior of m occurs at similar values of the ratio 
T/(co/a) ~ 0.6 for both S=l/2 and S=l models, above 
which a crossover to the T ^ J asymptotics m ~ 1/T 
occurs. The region where m becomes nonuniversal and 
no longer represents the spin wave mass is shaded 
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FIG. 5. m/T is plotted versus the ratio T/{co/a) for the 
S=l/2 (circles) and S=l (diamonds) single-plane models, as 
well as the model of two S=l/2 planes (triangles) coupled 
such that the system is at the critical point (as determined 
by Sandvik and Scalapino) and therefore should exhibit QC 
behavior at all temperatures. The difference between mod- 
els in the temperature range where scaling applies is evident: 
the S=l model is always renormalized classical, the two plane 
model is always quantum critical, and the S=l/2 model ex- 
hibits a crossover between the two around T = 0.5J. More- 
over, for S=l/2 and S=l models m/T deviates from the re- 
spective scaling temperature dependences at roughly the same 
value of T I (co I a) > 0.6 (shown as shaded area), as expected 
given the similar character of deviations from linearity in their 
spin wave spectrum (see text). For co, we use the spin wave 
theory results: co = 1.67 J for S=l/2, co = 3. 07 J for S=l 
[16], and co = 1.9 J for the two- plane model [17]. 

in Fig. 5. Expressed in terms of J, the range where uni- 
versal behavior of m(T) holds is below T ~ J for the 
S=l/2 model and below T ~ 1.9J for the S=l model. 

Now, we turn to the analysis of the spin wave velocity 
calculated in Sect. Ill, and plot c determined from the 
data at different temperatures using Eq.l3 in Fig. 6. 

Dashed lines show the known T = value of the spin 
wave velocity calculated using 1/5* expansions [16], 



: 2V2.JS(l + 0.1580/25 0.0216/(25)2) 

1.67J for S=l/2 
3.07J for S=l 



(23) 



which for S=l/2 compares favorably with the T=0 series 
expansion calculations [14]. 

For all temperatures studied the spin wave velocity 
turns out to be temperature dependent for both S=l/2 
and S=l models. Towards the lowest temperatures in 
the numerically accessible range, c approaches the ex- 
pected T=0 values given by Eq.(23). Since in the a- 
model the spin wave velocity is temperature independent 
except near the ultraviolet cutoff, which in our case is set 
by the scale of order c/a ~ 2\/2JS, the observed 
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FIG. 6. The spin wave velocity plotted as a function of tem- 
perature for the S=l/2 (a) and S=l (b) Heisenberg models. 
Dashed line shows the known T=0 values, co [18,16], and solid 
line is the result of 1/5 calculation by Kaganov and Chubukov 
[19], Eq.24. The temperature dependence for T > JS arises 
from nonuniversal lattice corrections. 



might also lead to a breakdown of scaling altogether by 
introducing deviations from linearity in the spin wave 
spectrum and appearance of spectral weight in modes 
not contained in the quantum nonlinear cr-model. We 
argue that this is not the case here. Indeed, from the 
fact that Eqs.(12,13) remain valid when c(T) is already 
temperature dependent, we conclude that the spin wave 
spectrum still remains linear even though c(T) deviates 
from its T=0 value. By the same token, we conclude that 
there are no modes with significant spectral weight other 
than those described by the cr-model. 

Another universal quantity which has not been dis- 
cussed in our paper so far is the overall prefactor, xQj 
in front of the scaling expression for x(q,'^) in Eq.(16). 
This prefactor, which is renormalized by short wave- 
length quantum fiuctuations, can be expressed in terms 
of No and ps for temperatures where corrections of or- 
der T/JS can be ignored. Where c is T-dependent, this 
quantity becomes T-dependent as well. However, quan- 
tum corrections to the prefactor due to l/S terms have 
not yet been calculated analytically. 

Now the reason why the quantum critical scaling fails 
to describe the temperature dependence of the correla- 
tion length is clear. The correlation length derived from 
Eq.(A3) is ^ = c/m. If one assumes that the spin wave 
velocity is temperature independent above T = 0.4 J, the 
quantum critical expression does not agree with the nu- 
merical data. Given the discussions above, we argue that 
one should consider not the temperature dependence of 
^, but that of m = c/^, because m is unaffected by the 
lattice corrections to c. Then, the crossover from RC 
behavior to QC behavior at temperatures around J/2 
becomes evident. 



V. CONCLUSION 



temperature dependent spin wave velocity clearly results 
from the quantum nature of spin, when the temperature 
becomes comparable to JS. 

The quantum corrections to the spin wave velocity 
originating from higher order l/S terms in the spin wave 
theory, were calculated by Kaganov and Chubukov [19]: 



c(T) = Co 



1 



CM 

8TrS 



T 
5J 



(24) 



where C(*) is the Riemann zeta function. The predictions 
of Eq.(24), evaluated for S=l/2 and S=l, are shown as 
solid lines in Fig. 6a, b, respectively. The agreement with 
numerical calculations is quite good. Despite higher or- 
der terms in the expansion in powers of T/JS, it appears 
that Eq.(24) describes the observed temperature depen- 
dence of c fairly well for both 5=1/2 and S = I. 

At this stage an important question is, does the tem- 
perature dependence of c destroy scaling altogether? One 
might expect that once lattice corrections become im- 
portant for quantities like the spin-wave velocity, they 



Magnetic properties of the doped superconducting 
cuprates may be closely related to the undoped parent 
compounds, such as La2Cu04 and Sr2Cu02Cl2, which 
are well described by the square lattice Heisenberg model. 
Motivated by this relationship, we have studied scaling 
and the role of lattice corrections in two dimensional 
square lattice Heisenberg antiferromagnets. 

In [3], Chakravarty, Halperin, and Nelson mapped 2D 
collinear quantum antiferromagnets onto the quantum 
nonlinear cr-model, and derived the corresponding phase 
diagram. When the T=0 ground state of the correspond- 
ing cr-model has Neel order, as is the case for the square 
lattice Heisenberg models, the low temperature scaling 
regime is classical (often called renormalized classical, or 
RC), where the characteristic energy scale ui <^T, while 
at higher temperatures a crossover to the quantum crit- 
ical (QC) regime may be observed, where w ~ T. The 
S=l/2 square lattice Heisenberg model, which describes 
La2Cu04, was shown to exhibit RC behavior at low tem- 
peratures [3]. As it has a substantial ordered moment. 
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it was earlier believed not to have a region of quantum- 
critical behavior at intermediate temperatures. 

More recently, Chubukov and Sachdev [5] argued that 
the crossover from the RC to nonuniversal microscopic 
behavior occurs through the intermediate QC regime at 
0.4J < T < 0.6J, with universal behavior extending up 
to T = 0.6J and nonuniversal behavior above that tem- 
perature. The importance of this at first sight minor 
difference concerning a fairly narrow temperature range, 
arises from the fact that upon doping the QC region 
was argued to expand rather rapidly, covering a much 
broader, and well accessible temperature range [5]. This 
scenario has important implications for our views of mag- 
netism in the superconducting cuprate materials [5,20]. 

In our earlier paper [12], we showed that the scaling 
functions for the equal time spin correlator, S(q), and 
the static susceptibility, x(ct)) in the S=l/2 Heisenberg 
model in the temperature range T = 0.6 J — J agree with 
those calculated for the cr-model in the QC regime. The 
same turned out to be true for the ratio 5'(q)/7"x(q). Our 
comparisons did not involve any adjustable parameters. 

In the present paper, we undertook a detailed study of 
these issues using the high temperature series expansions 
approach. The simultaneous study of S=l/2 and S=l 
models allows us to contrast the two cases and highlight 
the parameter-free agreement of the intermediate tem- 
perature behavior for the 5 = 1/2 case with the universal 
quantum critical predictions of the non-linear cr-model. 
For both 5=1/2 and S = I models, we find that the 
q-dependences of S(q) and x(ct) separately, as well as of 
their ratio, agree very well up to T ~ 0.6co/a with the 
(T-model predictions calculated separately using l/N ex- 
pansions (for l/N calculations, we utilized many of the 
results obtained in [5]). This allows us to calculate the 
temperature dependence of the spin wave mass, defined 
as m = c/^, for which both RC and QC predictions are 
known without adjustable parameters. 

For S=l/2, we find that m(T) is in good agreement 
with the RC prediction for T < 0.45J, while for 0.65J < 
T < J the results are very close to the QC prediction 
and far from the RC prediction, exhibiting a crossover be- 
tween these two regimes for 0.45 J < T < 0.65 J. Further- 
more, we find that for T > 0.5J, the lattice corrections 
cause a temperature dependence of the spin-wave veloc- 
ity c and affect the overall prefactor in front of the scaling 
expression for x(q,'^)- The new results strongly support 
our earlier conjecture [12] that lattice corrections are such 
that full q- and w-dependences remain universal in the 
QC regime for temperatures up to T ~ 0.6co/a ~ J: 

x(q,.)=x.<I>Qc(fc^,^), (25) 

where $'^'^(a;, y) is the scaling function of the QC regime, 
computable using, for instance, the approach of [5]. For 
the S=l model, where quantum fiuctuations are substan- 
tially less, we find that RC behavior directly crosses over 
to the high temperature interacting local-spin behavior 



without an intervening QC regime. Finally, the observed 
lattice corrections to c for both S=l/2 and S=l roughly 
agree in magnitude and temperature dependence with 
earlier calculations of quantum T/JS effects by Kaganov 
and Chubukov [19]. 

In their recent publication, Greven and coworkers [8] 
pointed out the agreement of the measured and numeri- 
cally calculated correlation length C(T) with RC predic- 
tion for all temperatures below T = J . In our work, we 
present an alternative theory which is not only consis- 
tent with the same data for C(T), but explains the q- 
and w-dependences of spin correlators as well. We show 
that detailed examination of the wavevector and temper- 
ature dependence of S(q) and xil) supports the picture 
of renormalized classical behavior for temperatures up to 
T = 0.45J, and a crossover to a quantum critical behav- 
ior around T = J/2. 

While the QC behavior for the spin wave mass occurs 
in our calculations in the temperature range T > 0.6J 
currently inaccessible to experiments in La2Cu04 and 
Sr2Cu02Cl2, many signatures of the QC behavior have 
been experimentally observed at lower temperatures. For 
instance, the nuclear relaxation rate, 1/Ti, in La2Cu04 
saturates above T = 600 A' ~ 0.4 J [7] at a value close 
in magnitude to the universal QC prediction [5]. This is 
consistent with the present study, as it has been shown 
[5,21] that for those quantities which in the RC regime 
depend on ^ through its logarithm only, such as the ra- 
tio T1T/T2G and the bulk susceptibility, xo, the RC-QC 
crossover should be shifted substantially towards lower 
temperatures, in agreement with the measurements of 
Xo [6] and TiT/T-zg [7] in La2Cu04. 

It has been proposed earlier that the temperature 
range of quantum critical behavior rapidly expands with 
doping, and that the scaling behavior described by the 
quantum nonlinear cr-model may be observed in doped 
antiferromagnets, and specifically in the high tempera- 
ture superconductors [5,20,12]. We hope that the study 
of the renormalized classical to quantum critical scaling 
crossover and of the origin of nonuniversal corrections to 
scaling in the insulator presented here, will be helpful in 
understanding magnetic behavior of the doped systems 
as well. 
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APPENDIX A: QUANTUM NONLINEAR 
(7-MODEL 

The mapping [3] of the low energy spectrum of the 
Heisenberg model to the quantum nonlinear sigma model 
has been extensively discussed in the literature and text- 
books (e.g., see [22]). Here we briefly discuss some as- 
pects of this mapping, and of the 1/A expansion in this 
model, as necessary for our purposes. We refer the reader 
to Refs. [3,5] for a far more detailed treatment. 

The action of the quantum nonlinear cr-model can be 
written as: 



S 



PsO 

2 



l/T 



dr 



1 



1, 



(Al) 



where the constraint = 1 captures the scattering ef- 
fects and gradient terms lead to the linear spin wave spec- 
trum. Here n is the antiferromagnetic order parameter, 
Pso is the bare spin stiffness and cq the bare spin wave ve- 
locity. The assumed ultraviolet cutoff is set by the lattice 
scale of the underlying Heisenberg model. 

This model has been extensively studied in recent years 
using renormalization group methods [3], Monte Carlo 
simulations [13], and 1/A expansions [5]. The results 
obtained using different techniques agree with each other 
quite well. In our calculations, we rely mostly on the 1/A 
expansions approach. In the next Subsection, we discuss 
the leading order of this expansion, namely, the N = 
cxD approximation, and then proceed with the subleading 
1/A corrections. 



1. N=oo Approximation 

The N = CO approximation is based on the replace- 
ment of the physical 0(3) quantum nonlinear cr-model, 
where the three-component order parameter n corre- 
sponds to three-component physical spin, by the O(cxd) 
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(T-model where the number of components of n is infinite. 
A consequence of such a replacement is the absence of 
damping. This makes the N = co approximation unsuit- 
able for such problems as the nuclear spin-lattice relax- 
ation or the neutron scattering, where correct description 
of damping is essential. Nevertheless, for our particular 
purposes of extracting scaling parameters from the data, 
already the N = co approximation is relatively accurate. 
It is further improved by taking into account the first 
subleading corrections, which are discussed in the next 
Subsection. 

The unit length constraint can be enforced using La- 
grange multiplier. A: 



S 



PsO 

2 



l/T 



dT 



;{drnf + {Vnf + i\ {n^ - l) 



(A2) 



The action is now quadratic in n which can therefore 
be integrated out, leaving a functional integral for the 
field A. In case of N = co, the A field does not fiuctuate 
around the saddle point value «(A), which in what follows 
is denoted as for the reason which becomes clear mo- 
mentarily. By performing Fourier transform of Eq.(A2), 
one obtains the N = co solution for the Matsubara fre- 
quencies [5]: 



A 



+ c^q^ -\- Lo'^ 



m 



(A3) 



where c is the spin wave velocity and it is now evident 
that m = c/^ has the meaning of the spin wave mass, i.e. 
the gap in the spin wave spectrum. 

The Matsubara correlator Xsil, ii^n), where ui„ = 
2'iTnT, is related to the physical real frequency response 
function by analytical continuation: uj iui„. It is easy 
to see that Eq.(A3) corresponds to a spectrum of un- 
damped magnons with a gap: 



A 



-\- c^q^ 



A 
2Z 



1 



1 



where the spin wave dispersion is: 



(A4) 



(A5) 



The imaginary part of x(?,'^), which describes dissipa- 
tion, is: 

(A6) 

Of course, for the physical 0(3) QNLcr model the damp- 
ing is finite and the peaks in x" around to = ifg acquire 



finite width. This effect is captured only in the terms 
beyond the N = 00 approximation. 

Our considerations for the S=l/2 and S=l Heisenberg 
models are based on the high temperature series data 
for the q-dependent static susceptibility of the staggered 
order parameter (n) field , Xs{l)j and equal time struc- 
ture factor of the same field, S{q). Using Eq.(A3), one 
obtains: 



x{(i) = xAq,i^n = 0) 



A 



(A7) 



+ 00 



Siq)=T J2 X.(9,*^n) = -tanh-l , ^^g^ 



n = — 00 



where 

A = Xs{0)xm\ = (m^ -Kc2g2)i/2 

It follows directly from Eq.(A7) that: 
X.(0) 



Xs{q) 



1+ — 9 for# = cxD. 



(A9) 



(AlO) 



Further, the ratio of Eq.(A7) and Eq.(A8) depends only 
on Cg/T: 



i^tanh-i f^l 
2T 12T\ 



Txsiq) 2T L2TJ 
We now introduce a function ')[x\ defined as: 
x/2 



tanh(a;/2) 



X. for X > 1. 



(All) 



(A12) 



which has the following asymptotic behaviors: ')[x'\ ~ 
2^3(a; - 1) for < a; - 1 < 1 and -^[x] ~ 2x for a; > 1. 
Then 



S{q) 



TXs(q) 



for^ = o«- (A13) 



2. 1/N Expansion 

In this Subsection, we discuss corrections to the N = 
cxD equations (A10,A13) for the physical case of # = 3. 
We utilize the results of the l/N expansion calcula- 
tions by Chubukov, Sachdev, and Ye [5] and find that 
Eqs.(A13,A10) remain quite accurate even for N = 3. 

We turn first to Eq.(A13). While the full q-dependence 
of l/N corrections have not been calculated so far, the 
results for small wavevectors as presented in [5] turn 
out to be quite sufficient for our purposes. In the 
asymptotic high temperature (quantum critical) regime, 
m/T is temperature independent, while Xs(0)/Xs(q) and 
S(q)/Txs(q) depend only on the product q^ ~ cq/T. 
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With straightforward application of the results of [5] , we 
obtain for small q: 



Xs{q) 



1 



Vl—] q 



TXs(q) 



U2 



T 



V2 



T 



(A14) 



(A15) 



where 



vi = l- 0.0042/#^ 0.9986, 
U2 = l + 0.1473/#^ 1.0491, 
V2 = l + 0.0382/#^ 1.0127 



(A16) 



(l/N terms are evaluated for N = 3). 

The l/N corrections to the coefficients vi, V2, which 
enter in parenthesis in front of in Eqs.(A14,A15), are 
essentially negligible. Higher terms in the expansion in 
powers of q"^ have not been calculated so far. However, 
once the leading terms have very small corrections, cor- 
rections to higher order terms in the expansion in q"^ 
should be of the same order of magnitude. The reason 
for this conclusion is as follows. Eqs.(A14,A15) can be 
considered as relationshps between the position of the 
closest to the origin singuliarity for imaginary wavevec- 
tors, (/sing = l/'f, and the expansion in powers of q"^ near 
q = 0. Were the higher order terms large, the position 
of the singuliarity expressed in terms of the coefficient 
in front of g^, namely, vi or V2, would change, which in 
its turn would lead to large corrections to vi,V2. Given 
the fact that such corrections are in fact small, and ig- 
noring the possibility of accidental cancellations for both 
Eqs.(A14,A15), we arrive to the conclusion that higher 
order q"^" terms are small and therefore Eqs.(A14,A15) 
are applicable for qa I, not only for the much narrower 
wavector range ^ 1. 

It is quite remarkable that the l/N correction to the 
coefficient U2 is small even when evaluated at # = 3. 
Indeed, the ratio S(Q)/Tx(Q) reflects the distribution 
of the spectral weight over frequencies. At the com- 
mensurate wavevector, the character of this distribution 
changes from a delta-function peak at flnite w = m to a 
smooth peak centered at w = as one goes from N = co 
to N = 3, which in general could have led to large l/N 
corrections to U2. The fact that the rigorously calcu- 
lated correction (A16) is in fact small shows that already 
the N = CO approximation accurately captures the re- 
lationship between the spin wave mass and the relative 
correlation strength at different Matsubara frequencies. 

Further, since the above corrections are essentially neg- 
ligible in the QC regime where they have been calcu- 
lated, it iseems reasonable to conjecture that they may 
not increase drastically as m/T decreases, in which case 
they should apply to the RC-QC crossover regime as 
well. The rigorous calculations to check this conjec- 
ture are prohibitively complicated and therefore have not 



been carried out. Instead, in the next Section we verify 
Eqs.(A14,A15) are in agreement with the numerical data 
for Heisenberg models for all temperatures studied, and 
thereby establish that the above conjecture is indeed cor- 
rect, which greatly simplifles futher considerations. 

To summarize, in this Appendix we derived surpris- 
ingly simple Eqs.(A14,A15,A16). Negligible (for vi, V2) 
or small (~ 5% for U2) l/N corrections to the N = co 
result indicate a well behaved large N expansion, and 
we expect the above equations to be accurate for the 
physical 0(3) model. The correspondence between the 
Heisenberg and sigma-models is such that Xsil,'^) for 
the (T-model is x(|ct — Ql,"^) for the Heisenberg model, 
where Q = (it, it). Taking this into account, we rewrite 
Eqs.(A14,A15) in the form appropriate for the Heisen- 
berg models, Eqs.(12,13), and use them in Section III for 
comparisons with the numerical data. 



APPENDIX B: CORRELATION LENGTH 

In this Appendix, we verify our resluts for m and c 
by comparing the ratio c/m with the independently cal- 
culated correlation length. We calculate the correlation 
length, deflned as the inverse rate of exponential decay of 
spin correlator at large distances, directly using our new 
"imaginary wavevector" method [23]. The idea of this 
calculation is based on an observation that when the pair 
correlation function decays exponentially up to power- 
law corrections, which is the case for two-dimensional 
Heisenberg models, the Fourier transform of spin corre- 
lator around the ordering wavevector, x(Q + ng(/), where 
iiq is a unit vector, should have a singuliarity on the 
imaginary q axis at (/sing = i/C when analytically contin- 
ued to complex q. The position of this singuliarity can 
be determined through high temperature series expan- 
sions using the methods described in [23]. In Fig. 7, our 
results for ^ for the S=l/2 Heisenberg model are plotted 
along with the earlier quantum Monte Carlo calculations 
by Makivic and Ding [10]; the agreement is excellent at 
all temperatures. 

Now we can verify the equality ^ = c/m, which holds 
by the deflnition of m. In Fig. 8, we plot our directly cal- 
culated ^ by a solid line, the ratio c/m determined solely 
from Eq.(12) as hollow circles, and the ratio of c and m 
with both determined seperately from Eq.(13) as solid 
circles. The agreement between all three sets of data is 
excellent for both S=l and S=l/2. This agreement not 
only indicates that our calculations of m and c are ac- 
curate, but also serves as another veriflcation that the 
(T-model corectly describes wavevector and frequency de- 
pendences of the dynamical susceptibility of the Heisen- 
berg models in the studied temperature range. 
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FIG. 7. The correlation length determined using our new 
method ol analytical continuation to imaginary wavevectors 
(solid line) plotted along with earlier Monte Carlo results or 
Makivic and Ding [10] (circles). 
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FIG. 8. The correlation length for the S=l Heisenberg 
model calculated directly using the analytical continuation to 
imaginary wavevectors (solid line), and independently using 
the mapping to the quantum nonlinear cr-model, from Eq.(12) 
(triangles) and Eq.(13) (circles). All three sets of data agree 
with each other for both S=l/2 and S=l. 
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